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Abstract 

We have calculated the resistivity of a GaAs/AlGaAs heteroj unction in the 
presence of both an in-plane magnetic field and a weak perpendicular com- 
ponent using a semiclassical Boltzmann transport theory. These calculations 
take into account fully the distortion of the Fermi contour which is induced 
by the parallel magnetic field. The scattering of electrons is assumed to be 
due to remote ionized impurities. A positive magnetoresistance is found as 
a function of the perpendicular component, in good qualitative agreement 
with experimental observations. The main source of this effect is the strong 
variation of the electronic scattering rate around the Fermi contour which 
is associated with the variation in the mean distance of the electronic states 
from the remote impurities. The magnitude of the positive magnetoresistance 
is strongly correlated with the residual acceptor impurity density in the GaAs 
layer. The carrier lifetime anisotropy also leads to an observable anisotropy in 
the resistivity with respect to the angle between the current and the direction 

of the in-plane magnetic field. 
PACS Numbers: 73.20.Dx, 73.40.Kp 
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I. INTRODUCTION 



Magnetotransport in a two-dimensional electron gas (2DEG) is an extremely rich and 
complex subject. In the more usual geometry with a magnetic field perpendicular to the 
plane of the 2DEG, the magnetoresistance exhibits Shubnikov-de Haas oscillations which 
are the low-field precursors of the quantum Hall effect. This behaviour is a consequence of 
the quantization of the in-plane orbital motion into Landau levels. The situation with a 
magnetic field in the plane of the 2DEG is quite different since the cyclotron motion now 
competes with the effects of the potential confining the electrons. This leads to an interesting 
modification of the electronic energy band structure which was recently analyzed in detail 
for typical GaAs/AlGaAs heterojunctionJlli. These calculations account for the observed 
depopulation of higher subbandsi which occurs with increasing parallel magnetic field. 

In comparison to the perpendicular field geometry, transport in the presence of a parallel 
field has received relatively little attention. One of the earliest studies is that of Englert 
et al.i who observed a pronounced positive transverse magnetoresistance when an in-plane 
magnetic field was tilted slightly out of the plane of the 2DEG. This interesting observation 
was not understood at the time. More recently, the same effect was studied in consider- 
ably more detail by Leadley et al.i for a number of GaAs/AlGaAs heterojunctions. The 
explanation of this effect is the primary focus of our work. 

A positive transverse magnetoresistance is known to occur in situations where conduc- 
tion is provided by different kinds of carriers, for example, a hetero junction with two or 
more occupied subbands. In accord with classical theoryi, samples with only a single oc- 
cupied subband do not show the effect (in fact, a negative magnetoresistance attributable 
to weak localization is usually observed). However the situation changes when the system 
is subjected to a parallel magnetic field which affects the electronic structure and leads to 
a field-induced anisotropy of the 2D Fermi contour. As noted by Leadley et al.i, this is 
analogous to the Fermi surface anisotropy found in metals and should similarly lead to a 
positive magnetoresistance^. Although the anisotropy of the Fermi contour is important, we 
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shall see that it is not the only effect contributing to the positive magnetoresistance. 

There have been a few theoretical considerations of this problem. The first calculation of 
transport in the presence of a parallel magnetic field was carried out by Tang and Butcheri'i 
who considered a model in which the 2DEG is confined by a harmonic potential. The 
electronic states and energy dispersion can be obtained analytically for this model which 
simplifies the solution of the transport problem. They further assume that the electrons 
are scattered by short range ^-function potentials and demonstrate within the framework 
of the Boltzmann transport theory that the conductivity of the electron gas is anisotropic 
with respect to the angle between the current and the parallel magnetic field. However they 
do not consider the additional effects of a perpendicular magnetic field. A similar model 
was used by Smrcka0 in an attempt to explain the observed positive magnetoresistance. 
However his calculation invokes a dc conductivity in the third dimension perpendicular to 
the plane of the 2DEG which has no obvious physical meaning. More recently, Stfeda et 
al.0 developed a semiclassical Boltzmann transport theory of the magnetoresistance within 
a constant relaxation time approximation. The theory assumes that currents flowing in a 
direction perpendicular to the in-plane magnetic fleld induce a transverse Hall fleld per- 
pendicular to the 2DEG as flrst suggested by Stern@. It is then argued that this field in 
turn leads to a correction to the in-plane conductivity which is the source of the positive 
magnetoresistance. However, it is unclear how this mechanism might account for the magne- 
toresistance when the current is flowing parallel to the in-plane fleld. In any event, since the 
argumentation is rather heuristic, a more rigorous explanation of the experimental observa- 
tions is needed. We have therefore performed a more thorough calculation of the transport 
properties in a tilted magnetic fleld, utilizing a realistic model of both the hetero junction 
electronic structure and the scattering of electrons from remote ionized impurities. We shall 
demonstrate that the inherent anisotropy of these properties is ultimately the source of the 
positive transverse magnetoresistance. 
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We perform our calculations for model heterojunctions which correspond to the G590 
series of samples studied by Leadley et al.l By modelling these particular samples we hope to 
minimize any differences in electronic structure which might otherwise obscure a comparison 
between theory and experiment. The design characteristics of these samples have been dis- 
cussed elsewherei. The important variables are the areal acceptor and donor densities {ua 
and no), which can be varied experimentally by means of the persistent photoconductivity 
effect. We shall deal specifically with two combinations, herein referred to as Sample 1 and 
Sample 2. Sample 1 has an acceptor density of ua = 1.03 x lO^^cm"^ and a donor density of 
no = 4.53 X lO^^cm"^, giving a total electron density of rix = 3.5 x lO^^cm^^. It corresponds 
to the unilluminated experimental sample referred to as G590a which has only one occupied 
subband. For Sample 2, the areal densities are ua = O.lx lO^^cm"^ and ud = 7.2x lO^^cm"^, 
resulting in a total electron density of tit = 7.1 x lO^^cm^^, corresponding to the experi- 
mental sample referred to as G590c. This sample is obtained from G590a by illumination 
and has two occupied subbands. A comparison of the theoretical and experimental zero-field 
mobilities provides one measure of how well we have been able to model these samples. We 
find 20.6m^/Vs and 37.9m^/V^s for Sample 1 and 2, respectively, which compare favourably 
with the experimental values of nAm'^ /Vs and 94.0m^/V^s. We are therefore confident that 
our model heterojunctions closely represent the actual experimental samples. 



II. THEORY 

A. Boltzmann Transport Theory 

In this section, we develop the Boltzmann transport equation for a 2DEG in the presence 
of a strong parallel magnetic field in the x-direction and a weak perpendicular component. 
The parallel field is taken into account quantum mechanically in the determination of the 
subband electronic structure (see Ref. H). As a function of the in-plane wavevector k, the 
dispersion of the subband energies is anisotropic and leads to noncircular constant 



4 



energy contours. An example of this behaviour is shown in Figure 1, where the contours 
are illustrated for Sample 2 for two parallel magnetic fields, one {B = 5T) for which two 
subbands are occupied and a second {B = lOT) for which only a single subband is occupied. 

The dynamics of the electrons in the presence of a perpendicular magnetic field compo- 
nent B_i will be described semiclassically. The transport properties arising from an external 
electric field E are then determined using the Boltzmann equation 



dt 



scat 



Ij E ■ Vk/n(k) - ^ Kk X B^) ■ Vk/„(k) (1) 



where /n(k) is the nonequilibrium distribution function of subband n and v„k = Vk-^nk/^ 
is the group velocity. Since the energy bands are neither parabolic nor isotropic in the 
presence of a parallel magnetic field, the electron velocity is not parallel to the wave vector 
k. Expressing /„(k) as 

/„(k) = /0(E„k)+^?n(k) (2) 

where /"(-Enk) is the equilibrium Fermi distribution and 5'ri(k) is the deviation from equi- 
librium, the linearized Boltzmann equation is given by 

5^n(k) 



icat 



dt 

The scattering term in (|^) is 

5^n(k) 



|j E ■ Vk/°(^nk) - ^ (v„k X B^) ■ Vk^7„(k) . (3) 



dt 



= ^[(7„,(k')-(7„(k)] Wr^A^M) (4) 

scat n'k' 



where if„„'(k, k') is the transition rate between the states nk and n'k'. It satisfies the 
detailed balance condition w„/„(k', k) = w„„/(k, k') which ensures the overall conservation 
of particle number. 

In the following we shall assume that the scattering is due to impurities. In this case, 
Fermi's golden rule gives the transition probability 

27r 



w„„,(k,k') = — |(nk|V(r)|n'k')|2 (5(E„k - ^n'k') (5) 



where V{r) is the screened impurity potentiaL The bar over the matrix element denotes an 
average over all possible impurity configurations. The important property for the present 
purposes is the energy- conserving delta function which corresponds to the elastic nature of 
the scattering. Making use of this dependence and defining the quantity (^n(k) by 

gM = ^^^a(k) (6) 

equation (H) can be reduced to 

T-\k) e„(k) = ^ Wnn'iK k') e„'(k') + eE ■ v„k + ^ (v.k X B^) • Vke„(k) (7) 

n'k' 

Here, r„(k) is the quantum lifetime 

r-^(k) = ^^„„,(k,k') (8) 

n'k' 

which represents the total unweighted probability of scattering from the initial state nk to 
all available final states, including states in other subbands. Due to the more complicated 
energy subband structure in the presence of a parallel magnetic field, w„„/(k, k') is not 
simply a function of k — k' and as a result, r„(k) is in general an anisotropic function of the 
wavevector k. 

The anisotropic energy band structure also complicates the solution of the Boltzmann 
equation. We find it useful to transform from the k-space variables to a curvilinear coordi- 
nate system defined by the energy variable e = Enk and an angular variable corresponding 
to the orientation of the vector k. To be specific, the position in k-space of a state in the 
nth subband is referred to the position of the minimum of the subband energy which is 
displaced from ky = because of the parallel magnetic field. This new origin is implicitly 
assumed in the following. The definition of these variables is illustrated in Figure 1. We 
have found that the use of the variable is more convenient in the present context than the 
phase angle 6 which naturally appears in discussions of cyclotron motion (Ref. Chap. 9). 

With this variable transformation, equation becomes 

r-i(e, 0) en(e, 0) = E /'^ 'P') ^nn'ie, 0, 0') ^^e, 0') 

ivr Jo 



+eE ■ v„(e, 0) - (^^) Vn{t, 0)t„ ■ Vk^nle, 4>) (9) 

where t„ is a unit vector pointing in the direction of cyclotron motion on the constant energy 
contour of the nth subband at the point k. The angular transition rate at energy e appearing 
in d^) is defined by 



Pnn'{En^,<f),<P') = Aj^ |(nk|\/(r)|r2'k')|2 (10) 



in terms of which the quantum lifetime is given by 
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V'(e,0) = E:^ f^d^'Jr.ie,^')Pnn'{e,^,<P'). (11) 
The integrals in (D and (pJ]) contain the quantity 

■^"''■^'' .n-lv'g kl '''' 
m I Vk-c^nk • k| 

which is a dimensionless form of the Jacobian of the variable transformation. In this form it 
is just the ratio of the free particle velocity hk/m* to the component of the actual velocity 
v„(k) in the direction k. The deviation of this function from unity reflects the asymmetry 
of the constant energy contours in the presence of a parallel magnetic field. 

The solution of (P) is easily obtained using a Fourier expansion in the angular variable, 

i.e. 

UeA) = T.^t\^)^'""^- (13) 

m 

Multiplying (|^) by e"*"''^ J„(e, 0) and integrating over yields the set of linear equations 

E - i?i:r"^'^ + tmu'j^^^S^n] ai?'^ = (14) 



m'n' 



with 



4"^^ = ^/^rf0e— '^J„(0)r-^(0), (15) 
/TT JO It: Jo 



and 

bt^ = ^ r rf0e-*'"^^(0)eE ■ v„(0) (17) 
Ztc Jo 

In these equations, the dependence on the energy variable e has been suppressed. We note 
that the free electron cyclotron frequency lu^ = appears in equation (0), and not the 
actual cyclotron frequency defined by uj~^ = § where the line integral extends over 
a constant energy contour. 

In the absence of a parallel magnetic field the energy dispersion is isotropic, Jn{4>) reduces 
to unity and the angular transition rate is only a function of the difference — </)'. As a 
result, T~^(0) becomes independent of (p, b'^^ is proportional to 5m,±i and a^^^^ are the 
only non- vanishing expansion coefficients. In this limit, one recovers the usual form of the 
multisubband transport equations0. 

In the present situation however, the equations defining the Fourier expansion coefficients 
are coupled. In practice, the expansion is truncated at some finite number of terms and the 
set of equations in is solved with the inhomogeneous term taking on two possible values, 
one corresponding to an electric field in the x-direction and one in the y-direction. Once 
these solutions are known, the transport current density is obtained from 
2e 

J = -^^^n(k)v„k 

nk 

Since the Fourier coefficient a^™^ depends linearly on the electric field, we can define the 
conductivity tensor as 

- s s £ * (-^) ^ r ^^•'"<^- ^>""''<'- ^''*"* • 'i*" 

In the B± limit, one can show that a^^ is diagonal, but (^yy As a result of the 

distortion of the Fermi contour, the current flow is not parallel to the applied electric fleld 
except for the special cases when the electric fleld is either parallel or perpendicular to the 
in-plane magnetic fleld. 
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We finally comment on whether or not the transverse Hall field proposed by Stfeda et ahllil 
is a possible mechanism for the positive magnetoresistance. The nonequilibrium distribution 
function in (|^) implies a spatial redistribution of the electronic charge in a direction normal 
to the interface given byi 

5n{z) = 2Y,g^{y)\(l^^,^{z)\' (20) 

rak 

where (pnkyiz) is the subband wavefunction in the presence of the in-plane magnetic field. 
This redistribution gives rise to a transverse Hall field and a Hall potential which is estimated 
to be of the order of microvolts under typical current carrying conditions.0'i These potentials 
are small on the energy scale of the confining potentiafl and will lead to small changes in 
the energy band dispersion. Corrections to the Boltzmann equation (for example, as a result 
of changes in the electron velocity) are therefore of higher order in the applied electric field 
and, at the level of the linearized Boltzmann theory, can be safely neglected. A transverse 
Hall field is therefore not relevant to in-plane transport. 

B. Evaluation of the Scattering Elements P„„/(k, k') 

We now turn our attention to the calculation of Pnn' (k, k') in (|1^) for a situation in 
which scattering is due to remote ionized impurities. Using the form of the subband states 
in the presence of a parallel magnetic field, the potential matrix element has the form 

(nk|V(r)|n'k') = jj rf#„,^(z)0„,fc, (z)\/(z, q) (21) 

where q = k — k' and 

V{z, ci) = j dse-'^^-^Vis, z) (22) 



is the 2D Fourier transform of the scattering potential. We note that the evaluation of ( |2TD 
for the parallel field configuration is considerably more involved because of the fc^-dependence 
of the subband states. As a result, the potential matrix elements do not simply depend on 
the momentum transfer q which appears in the Fourier transform of the impurity potential. 



We assume that the impurity potential arises from a distribution of remote ionized impu- 
rities located at sites (Rj, Z^). These impurities (of charge e) give rise to the bare electrostatic 
potential 



where a two dimensional Fourier transform has been taken in the plane of the 2DEG and 
we have assumed that the position of the charged impurities is such that Zi < z. k is 
the dielectric constant of the material in which the 2DEG is imbedded. The quantity /(q) 
determines the amphtude of the exponentially decaying Fourier transform and contains all 
information regarding the spatial distribution of the ionized donor impurities. 

The external impurity potential is screened by the 2DEG and it is the final screened po- 
tential which is responsible for the electron scattering. We shall account for the screening at 
the level of the random phase approximation in which the electrons respond self-consistently 
to the ionized impurity potential. Defining ^^(z, z', q) to be the 2D Fourier transform of the 
independent particle density response function, the change in electron density due to the 
impurities is given by 



where 4>^°*{z, q) is the total electrostatic potential acting on the electrons. This potential is 
given by 




(23) 




(24) 



(/)*°*(z,q) = 0-*(z,q) + 0-'^(^,q) 



(25) 



where the second term is the induced electrostatic potential 




(26) 



arising from the induced electron charge density. 
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In the absence of a parallel magnetic field, the subband envelope functions are It- 
independent and the calculation of the density response function is relatively straightforward. 
It takes the form 



X%z,z',cD = j:Mz)Mz')xi{q) (27) 

where the index (3 represents a pair of subband indices {n,n') and ff^{z) is the product 
(pn{z)(j)n'{z) of two subbaud envelope functions. The factor is just the intersubband 

response function for an ideal 2DEG which depends on the subband structure through the 
subband energy levels and the position of the Fermi level. Analytic expressions for x^^Iq) 
are available (see, for example, Ref. plSf ). 



Once a parallel magnetic field is present, the subband states acquire a non-trivial ky de- 
pendence and x^i^j ^' 1 q) cannot be obtained analytically. To avoid an excessive and largely 
unnecessary amount of numerical work, we shall make a simple approximation motivated 
by the following observation. Although the individual subband states are strongly modified 
by the parallel field, depending on the value of ky^ the overall electron density distribution 
is relatively insensitive to the field. This is illustrated for Sample 2 in Figure 2, where we 
compare the ground state density for B\\ = OT to the density at the relatively high field of 

= lOT. As can be seen, the effect of the field is minor, even for this case in which the 
second occupied subband is depopulated by the magnetic field. One would therefore expect 
the screening of the impurity potential to take place in the presence of the parallel field in 
much the same way as in the zero- field limit. We therefore adopt the physically reasonable 
approximation of screening the impurities by the zero- field response function given by (pT]). 

However, one complication must be addressed: higher occupied subbands will depop- 
ulate with increasing parallel magnetic field. Since different subbands screen differently, 
maintaining the zero field populations in the calculation of x^{z,z',q) introduces an er- 
ror which can be avoided in the following way. The subband populations ni{B) are first 
determined from a fully self-consistent electronic structure calculation which is known to 
reproduce the observed field dependence quite accura tely.0 We then make use of this in- 
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formation in the calculation of x${(l) by simply defining fictitious subband energies Ef to 
ensure that the subband densities are given correctly, that is, 

27mi{B\\) = Ep - . (28) 

This is sufficient to define the intersubband response function for the occupied subbands. In 
the cases we have dealt with, a maximum of two subbands are occupied and we have therefore 
truncated x${(l) to a two-by-two matrix corresponding to the lowest two subbands. Once 
the second subband is depopulated, the second subband still contributes to the screening by 
providing final states to which the first subband electrons can be excited. In this situation 
the second subband is positioned relative to the Fermi level according to the self-consistent 
calculation, while the first subband energy is still determined by (pSj). In practice we have 
found that these refinements have only a slight effect on the calculated scattering matrix 
elements since the lowest subband holds most of the electrons and contributes most of the 
screening. Nevertheless, since relatively little additional effort is required to determine these 
corrections, we have retained them in all of the screening calculations. 

With this prescription for x${(l) and z',q), the screened impurity potential is ob- 
tained from the set of equations (|23D-(|26D. Expressing (q, z) as 

0*''*(^,q)= (^^^ I{q)J{z,q), (29) 

we have 

J{z,q) = e""^ - (^^^ Y.xiiq) J dz' j dz" e~'^\^~^'\ fp{z')fp{z")J{z\q) . (30) 

We note that this integral equation has a separable kernel and it can therefore be reduced 
to a matrix problem. Multiplying (|30| ) by fa{z) and integrating over z gives 

Uq) = [ dzUz)e-'^^ -Y.Fo.Mxl{q)JM (31) 

where 

J,(g) = j dzUz)Jiz,q) (32) 
12 



Pnn'iK^') = ^1 ( J dz<P^ky{z)(pn'k'^iz)J{z,q)\ -|/(q)|2 (34) 



and 

F.p{q) = (^^) JdzJ dz'f^{z)fp{z')e-'^\^-^'\ (33) 

is a Coulomb form factor. As mentioned earlier, we retain two subbands in our calculations, 
so that the screening effect of the second subband is still present even when it is depopulated. 

The solution of ( |3TD is substituted into (^) to obtain the 2;-dependence of the screened 
impurity potential required in the calculation of the scattering matrix elements. In partic- 
ular, the angular transition rate takes the form 

m* I ( 27re^\ /",, , ^ ,|2l 

nq 

where the bar over the final factor denotes a configuration average. For uncorrelated impurity 
positions, we have 

^WF = Jdzp,iz)e'^^ (35) 

where pi{z) is the average spatial density of ionized impurities. For the special case of a 
delta-doped layer with Pi{z) = n^S^z + s), (^) becomes n^e"^''*, showing that the effect of 
the impurities diminishes exponentially with the setback distance s. 

We can compare this result with the idealized situation of short range impurity scattering 
considered by Tang and Butcheri. For a sheet of 5-function scatterers located in the plane 
z = zq, we find P„„/(k, k') oc |0nfcj;(-2o)Pl0n'fc;(-2o)P) which must be integrated over zq to 
correspond to a uniform distribution of scatterers throughout the region of the 2DEG. This 
scattering rate is independent of kx and k'^ which, as noted by Tang and Butcheri, simplifies 
the solution of the Boltzmann equation in the case of E || By. For the more realistic situation 
of remote charged impurities this simplification does not arise and the E || B|| case requires 
a treatment similar to that of the E ± Bn case. 



III. RESULTS AND DISCUSSION 



All of the interesting magnetotransport effects in the presence of an in-plane magnetic 
field ultimately arise from the distortion of the Fermi contours illustrated in Figure 1 and 
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the associated behaviour of the subband wave functions. That this could lead to a posi- 
tive transverse magnetoresistance was already appreciated by Leadley et al. who viewed 
the electrons on the distorted Fermi contour as different kinds of carriers having different 
mobilities. The usual treatment of parallel transport^ would then give rise to a positive 
magnetoresistance. They attributed variations in the mobility to a k-dependent effective 
mass which is due to the altered band structure. However, within a Boltzmann transport 
theory it is the electron velocity which emerges as the important dynamical variable, and 
the effective mass appears only when the energy dispersion is strictly parabolic, which is 
not the case in the present situation. In any event, a much more significant factor is the 
variation of the electronic scattering rates around the Fermi contour. In the following we 
shall try to differentiate between the kinematic effects which arise from the energy band 
structure and the dynamical effects associated with impurity scattering. 

Figure 1 shows an example of the Fermi contour anisotropy for a sample in which two 
subbands are occupied at zero field. As can be seen, the distortion of the Fermi contour of 
the first subband increases with parallel field strength and eventually takes on the shape of 
an egg. As this distortion is developing, the energy separation between the first and second 
subbands increases and leads to the depopulation of the second subband. (For Sample 2, 
complete depopulation occurs at approximately B\\ ~ 5.6T.) Interestingly, the minimum in 
the second subband also displaces relative to that in the first, so that the pocket of second 
subband states approaches the first subband Fermi contour with increasing field. This too 
will be seen to have an important effect. However, we shall begin by considering the simpler 
situation in which only a single subband is occupied, either because the electron density is 
low or because the field is sufficiently high to have depopulated the higher subband. Since 
our calculations are done for zero temperature, only the states at the Fermi energy are 
relevant and we can restrict the solution of (|^) to e = Ep. It should be understood that all 
quantities are calculated at this energy. 

One measure of the Fermi contour anisotropy is the Jacobian defined in (P^. In Figure 3 
we plot the inverse of this quantity as a function of the angular position around the Fermi 
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contour, together with the ratio of the magnitude of the velocity to the free-electron velocity 
hk/m*. We recall that the k- vector in this context is defined with respect to the position 
of the subband minimum and not the more usual k-space origin. One consequence of this 
definition is that there is only a small difference between the curves in Figure 3, indicating 
that the normal to the Fermi contour does not deviate much from the direction of k. In this 
respect, the Fermi contour is still rather circular in nature, despite its visual appearance. 
On the other hand, the four-fold variation of Jf^(0) around the contour is showing that the 
modified energy dispersion is having a dramatic effect on the electron velocity. There are 
two effects coming into play. The first is the overall elongation along the fcy-axis which is a 
result of a magnetic field enhancement of the effective mass in this direction. As discussed 
previously^, this leads to an enhancement of the electronic density of states above the ideal 
2D value. It is this effect that accounts for the decrease in J{^((f)) near = tt. However this 
same effect is swamped near = by a much larger reduction of the velocity coming from 
the flattening of the energy bands. Only at = 7r/2 is there no effect of the Fermi contour 
anisotropy, as the dependence of the energy on is unchanged by the parallel magnetic 
field. 

The difference in behaviour of the energy dispersion near = and = vr is a reflection 
of the fcy-dependence of the electron wavef unctions. Electrons with ky negative (0 tt) 
experience a magnetic potential which drives them into the interface, whereas states with 
positive ky (0 — > 0) are pushed away from the interface. The latter is a stronger effect, since 
the confining potential in a direction away from the interface is much softer, particularly 
when the background acceptor density is low. These differences are illustrated in Figure 4 
which shows the centroid of the subband probability densities as a function of ky. One can 
see that there is only a slight change in the centroid position as the wavefunctions are driven 
into the interface (7r/2 < < tt), but there is a much larger variation for those states being 
pushed in the opposite direction (0 < < 7r/2). 

The change in mean position has a dramatic effect on the angular scattering rate (P^. 
It is particularly sensitive to the wavefunction position since the bare scattering potential is 
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an exponential function of the distance from the impurity layer. Furthermore, those states 
which are further from the interface experience the full screening effect of the 2DEG which 
is interposed between them and the ionized donor impurities. The angular scattering rate 
depends on the two angular variables and 0' which specify the orientation of the initial 
and final wavevectors, respectively. Part of this angular dependence arises from the impurity 
form factor |/(q)p which depends on the momentum transfer q = k' — k. This factor is 
common to the zero field limit and we therefore choose to plot a normalized transition 
rate which has this factor removed in order to isolate the new effects associated with the 
dependence of the scattering matrix elements on the subband states. In particular, for the 
case of a single subband we consider the quantity 



where A0 is the angle through which the electron starting at the point is scattered. 
Because of the anisotropy of the Fermi contour, Pui<P, 4> + is not a symmetric function 
of the scattering angle, except at the special points = and = tt. More generally, we 
have the symmetry Pii(0, 0') = -Pii(27r — 0, 27r — 0'). 

Figure 5 shows the normalized scattering rate as a function of the starting angle for 
several scattering angles. For a given scattering angle A0, the momentum transfer q is 
approximately constant and the dependence on is an indication of the strong anisotropy 
arising mainly from the fcj^-dependence of the subband wavef unctions. The curve for A0 = 
corresponds to the limit of small angle scattering and shows a strong maximum at = tt. 
At this angle, the subband states are closest to the impurities where the screened potential 
J{z,q) is relatively large. As — > (or 27r), the subband states move away from the 
impurities and the scattering rate diminishes accordingly. A similar behaviour is seen for 
the other scattering angles, although now an asymmetry with respect to ±A0 is evident. 
The decreasing magnitude of the scattering rate with increasing scattering angle is partly 
due to the momentum transfer dependence of the screened potential J{z, q) which decreases 
with increasing q. 



Ai(0,A0) 



A 



Pn(0,0 + A0). 



(36) 
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It is also of interest to consider the anisotropy of the quantum hfetime defined in (|]). 
Given that the scattering rate drops off rapidly with increasing momentum transfer q, the 
inverse quantum lifetime is dominated by small angle scattering. In Figure 6 the anisotropy 
of the quantum lifetime is illustrated for Sample 1 at a field of lOT. The lifetime shows a 
four-fold variation between its maximum at = and its minimum at = vr, the latter 
occurring when the angular scattering rate has its maximum, as shown in Figure 5. To isolate 
the effects of the shape of the Fermi contour itself, we can perform a model calculation in 
which the actual subband wavefunctions are replaced by their zero-field limit. The lifetime 
anisotropy in this case is shown by the dashed line in Figure 6. Since the angular scattering 
rate in (|34D is now only a function of the momentum transfer, its anisotropy is relatively 
weak and the lifetime anisotropy is dominated by the Jacobian factor in ([TTD . As a result, 
the dashed curve in Figure 6 mimics the behaviour of Jf ^(0) in Figure 3. It is clear from a 
comparison of the two curves in Figure 6 that the full lifetime anisotropy is coming mainly 
from the magnetic field dependence of the subband wavefunctions. The variation with 
parallel magnetic field strength is shown in Figure 7. The lifetime at zero field is of course 
isotropic and the anisotropy about this value is seen to increase approximately in proportion 
to the parallel field. 

To illustrate the effect of multiple subband occupancy, we have performed similar life- 
time calculations for Sample 2 at an intermediate field (5T) where both subbands are still 
occupied. In this instance, we can define multi-band quantum hfetimes by the equation 

r-^' (0) = r-^' (k) = E ^nn' (k, k') . (37) 

k' 

These lifetimes represent the unweighted probability of an electron scattering from the state 
k(0) in subband n to all other states in subband n' . The two-band quantum lifetimes 
presented in Figure 8 can be seen to exhibit quite different anisotropies. For rn, we basically 
have the same behaviour shown in Figure 7 for the single-subband sample, but the magnitude 
of the anisotropy is significantly larger because of the smaller acceptor density in this sample. 
T22 has a similar behaviour, although the anisotropy is reduced because of the smaller range 
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of ky values in the second subband. For T21 we see a slight increase in the lifetime as — > vr. 
This arises since the second subband states near = vr tend to experience a slightly larger 
average momentum transfer when scattering into the first subband than do the states near 
= 0. However, the anisotropy of this intersubband lifetime is again small due to the small 
size of the pocket of second subband states. The variation of T12 is much more dramatic, 
covering almost three orders of magnitude. This intersubband lifetime shows a pronounced 
minimum near ~ vr/S which is the region where the first Fermi contour comes closest to 
the second subband. The reduced average momentum transfer for these states results in a 
scattering 'hot-spot' at which the lifetime is relatively small. 

In the zero field limit, the single subband transport lifetime differs from the quantum 
lifetime by the appearance of an additional (1 — cos0) weighting factor in the integrand of 
(P). In the presence of a parallel magnetic field, there unfortunately is no similar expression 
which can be used to define a transport lifetime. Instead, one must deal directly with 
the non-equilibrium distribution function as determined by the solution of the Boltzmann 
equation. This distribution function can always be expressed in the form 

e„(0) = E ■ A„(0) (38) 

which introduces the vector mean- free-path A„(0). In general, the mean-free-path is not 
parallel to the velocity vector at the point and for this reason the conventional definition 
of a transport lifetime is not appropriate. Nevertheless, for a given direction of the electric 
field, it is possible to parameterize the mean-free-path in terms of a lifetime. For example, 
for E = i?x (i.e., parallel to the in-plane magnetic field) we can write 

AnM ^ VnMrnM) ■ (39) 

Since the zeroes of ^n(0) coincide by symmetry with those of Vnx{4>) when the electric field is 
in the x-direction, r*^(0) is a well-defined quantity. This is not the case if the electric field 
is oriented in the y-direction and a similarly defined lifetime would exhibit singularities at 
certain points on the Fermi contour which have no physical significance. We shall therefore 
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simply use t*^(0) as a convenient parameterization in order to visualize the anisotropy of 
the solution to the Boltzmann equation. 

Figure 9 shows the transport lifetime of Sample 1 defined according to (p9D, normalized by 
the zero field transport lifetime, at three magnetic fields. For lower fields we see a monotonic 
angular variation similar to that of the quantum lifetime, with those states at — having 
a 50 percent greater lifetime than those at — vr. The states near = 7r/2 which contribute 
the most to the current have a lifetime very similar to the zero field lifetime. As the field 
strength increases, the transport anisotropy changes in character, with additional structure 
appearing near = 0. More importantly, the average lifetime decreases with increasing 
field and results in a positive parallel field magnetoresistance. This effect is illustrated in 
Figure 10 which shows the magnetoresistance for the current parallel and perpendicular 
to the in-plane magnetic field. The magnetoresistance is approximately parabolic and of 
a sizable magnitude for field strengths of the order of lOT. In addition, we see that the 
magnetoresistance is larger for J ± By than for J || By, in qualitative agreement with 
experiment^, although it should be noted that the experimental anisotropy seems to vary 
from sample to sample in an unpredictable way. 

We finally consider the transverse magnetoresistance which arises with a perpendicular 
component of the magnetic field. This component induces a cyclotron motion of the electrons 
around the Fermi contour and leads to an averaging of the anisotropic scattering rates. 
Figure 11 illustrates the behaviour found for a series of parallel magnetic field strengths for 
a situation in which two subbands are initially occupied. The model parameters in this case 
are ua = 0.3 x lO^^cm"^ and n^j = 7.4 x lO^^cm"^ {riT is still 7.1 x lO^^cm"^) instead of the 
values n^i = 0.1 X lO^^cm"^ and no = 7.2 x lO^^cm"^ used for Sample 2. This adjustment 
was made since the two subband mobilities for Sample 2 are almost identical and no positive 
magnetoresistance is found for zero parallel field. By increasing the acceptor concentration 
the calculated magnetoresistance correponds more closely to that observed experimentally 
and a more meaningful comparison of the parallel field dependence can be made. In addition, 
with the adjusted parameters the second subband depopulates at approximately 4T, which 
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is in better agreement with the G590c sample of Leadley et aLH 

The -Bji = OT curve in Figure 11 is the usual positive magnetoresistance associated with 
the different carriers in the two subbands (although it should be stressed that the carriers are 
not independent since intersubband scattering is included^) . As the parallel magnetic field 
is increased, both the depth of the magnetoresistance dip, and the field at which saturation 
occurs, increase. This is in good qualitative agreement with the observations of Leadley et 
al.i on their G590c sample, although the saturation field we find is approximately twice the 
value observed. Since saturation occurs when uJcT ~ 1, the discrepancy can be explained in 
terms of the difference between our calculated mobility of fi = 37.7 m'^/Vs and the experi- 
mental value of 94.0 m'^/Vs. Once the depopulation field of about 4T is exceeded, we find 
that the saturation field decreases suddenly to a value below the zero parallel field value due 
to the elimination of intersubband scattering and then remains constant, in marked contrast 
to the behaviour found below the depopulation field. The increase of the magnetoresistance 
with increasing parallel field in the single-subband limit is simply a consequence of the in- 
creasing lifetime anisotropy shown in Figure 7. This behaviour is in excellent qualitative 
agreement with that observed. 

In Figure 12 we give an example of the kind of anisotropy that could be expected in 
the perpendicular field magnetoresistance. Figure (a) is for the same sample discussed in 
Figure 10 and shows a magnetoresistance which is only slightly larger for J || By than for 
J _L B||. Figure (b), corresponding to the -By = 5T curve in Figure 11, shows a much 
larger anisotropy with the magnetoresistance still being larger for J || By than for J ± 
B||. This relative magnitude of the two magnetoresistances persists over a range of sample 
parameters and is opposite to the relative magnitude oberved in one particular sample.0 
We have no explanantion for this difference. It should also be noted that the experimental 
magnetoresistance tends not to saturate, but passes through a maximum before falling at 
higher magnetic fields.i'0 

To further quantify our results we show in Figure 13 the normalized magnetoresistance 
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for a range of acceptor densities with the total electron density fixed at nj- = 7.1 x lO^^cm"^. 
A decrease in ua can be induced by illumination. In agreement with observation, the 
magnetoresistance increases with decreasing ua] Figure 14(a) shows this dependence in a 
different way. This variation reflects the effect of illumination on the confining potential. 
As Ha decreases, the confining potential becomes softer and the subband states are more 
strongly affected by a parallel magnetic field. In fact, we have found that Ap^x/ Po is far more 
sensitive to variations in ua than in any of the other material parameters. For example, 
Figure 14(b) illustrates the effects of varying the total electron density tit while keeping 
the parallel magnetic field constant. For low values of ha, an increase in the total electron 
density translates into a reduction of the anisotropy in the system and a lower Ap^x/po- 
This is mainly due to an increase in the confinement of the self-consistent heterostructure 
potential with increasing rij-. For higher values of ua, the acceptors are themselves providing 
a stronger confining potential, and the effect of varying tit is greatly reduced. 

A comparison of our results in Figure 13 with the experimental data of Leadley et al. 
permits an estimate of the acceptor densities in the experimental samples. We find that ha — 
0.35 X 10^^ cm~^ is appropriate for the illuminated sample (G590c), while ua — 1-2 — 1.4 x 
10^^ cm~^ for the dark sample (G590a). These values indicate that the change in acceptor 
density as a result of illumination is Aha — 1 x 10^^ cm~^, which is consistent with the 
previously determined valuei. Interestingly, our estimates suggest that a significant residual 
acceptor density remains after illumination, in contrast to what is sometimes assumed. If 
the acceptor density were eliminated completely, our calculations would yield a much larger 
magnetoresistance than observed. 



IV. CONCLUSIONS 



We have performed extensive calculations of the transport properties of a 2DEG in the 
presence of parallel and perpendicular magnetic field components. The in-plane component 
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leads to a significant perturbation of the electronic structure which manifests itself as a 
distortion of the Fermi contours and a magnetic field depopulation of higher lying subbands. 
These features are one source of the anisotropy that emerges in both the quantum and 
transport hfetimes. We have also analyzed the impact of the parallel magnetic field on 
the scattering of electrons from the long-range potential due to remote ionized donors. The 
scattering at different points on the Fermi contour is found to be highly anisotropic as a result 
of the field dependence of the subband wavef unctions. This detailed scattering information is 
then used in a Boltzmann transport theory which treats the perpendicular field component 
semiclassically. Our general solution of the Boltzmann transport equation includes fully 
the effects of the electronic structure, intersubband scattering and the anisotropy of the 
nonequilibrium distribution function. Application of the theory to a realistic hetero junction 
is found to yield results which are in good qualitative agreement with experiment. In 
particular, our calculations account for the observed positive transverse magnetoresistance 
and its detailed dependence on the magnitude of the in-plane magnetic field in both the 
one- and two-subband regimes. It is clear from our calculations that a careful treatment 
of the electronic structure, long-range impurity scattering and transport behaviour are all 
needed to obtain a complete understanding of the experimental results. 
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FIGURES 

FIG. 1. Constant Fermi energy contours for Sample 2 with a parallel magnetic field of (a) 5r 
(two-band occupancy) and (b) lOT (single-band occupancy), kx and ky are in units of inverse Bohr 
radii ((a*)~^). The angular variable <p defines the position on the Fermi contour with respect to 
the position of the first subband minimum, fcj"*"; a similar angular coordinate is defined for the 
second subband, but is not shown. 

FIG. 2. The spatial density profile of the 2DEG in Sample 2 as a function of distance from the 
interface, for zero magnetic field (solid line) and for a parallel magnetic field of lOT (dashed line). 

FIG. 3. The inverse Jacobian m*v„ • 1^/hk, along the Fermi contour of the first subband, for 
Sample 2 in the presence of a lOT parallel magnetic field (solid line) . Also included for comparison 
is the function m*Vn/^k (dashed line). 

FIG. 4. The ccntroid of the first subband wave function relative to the interface, as a function 
of angular position on the Fermi contour. The results are for Sample 2 in the presence of a lOT 
parallel magnetic field. 

FIG. 5. The angular scattering rate P\\{(j), /S.<j)) of the first subband, for Sample 2 with a lOT 
parallel magnetic field. The scattering rate is normalized by the small-angle scattering rate at 
(f) = TT. The curves are labelled by the scattering angle A^, with the solid curves corresponding to 
positive scattering angles and the dashed curves to negative angles. 

FIG. 6. The first subband quantum lifetime, as a function of the angular position on the Fermi 
contour, for Sample 1 in the presence of a lOT parallel magnetic field. The solid line is the full 
calculation using the fe^-dependent wavefunctions, while the dashed line has been computed using 
only the constant zero field wavefunctions as described in the text. 

FIG. 7. As in Figure 6, but for a range of parallel magnetic field strengths. 
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FIG. 8. The multiple-subband quantum lifetimes Tnn'{4>) foi" Sample 2 in the presence of a 5T 
parallel magnetic field, as a function of the angular coordinate. 

FIG. 9. The effective transport lifetime of Sample 1 as a function of the angular coordinate, for 
an applied electric field parallel to the in-plane magnetic field. The lifetime has been normalized 
by the zero-field transport lifetime of the first subband. The curves are labelled by the parallel 
magnetic field. 

FIG. 10. Variation of the magnetoresistance as a function of the in-plane magnetic field. The 
labels XX and yy correspond to J || By and J ± By, respectively. 

FIG. 11. The transverse magnetoresistance for a series of parallel magnetic field strengths. The 
parameters characterizing this sample are described in the text. The second subband depopulates 
at a parallel magnetic field of approximately AT. The curves have been offset arbitrarily for clarity, 
but the scale of resistance variation is the same for each. 

FIG. 12. The transverse magnetoresistance for J || By {xx) and J _L By {yy). (a) and (b) 
correspond to Samples 1 and 2, respectively, both at a field i?y = ST. 

FIG. 13. The normalized transverse magnetoresistance as a function of parallel magnetic field 
strength. The different curves are labelled by the acceptor density n^, in units of 10^^ cm~^. All 
samples have the same electron density of 7.1 x 10^^ cm~^. The low-field termination of the curves 
occurs at the point where the second subband becomes populated. 

FIG. 14. (a) The normalized magnetoresistance as a function of ua at three different parallel 
magnetic fields, for the same set of samples as in Figure 11. (b) As in (a), but for parallel magnetic 
field of lOr and varying electron gas density, nr, in units of 10^^ cm~^. 
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